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Предложена модель соединения материалов с использованием синтеза в твердой фазе и проведена оценка напряжений, воз- 
никающих в зоне реакции в процессе соединения материалов в условиях однородного прогрева рассматриваемого образца или 
поджигания с торца. Результаты численного расчета показали, что теплофизические свойства соединяемых материалов и соеди- 
нительного состава существенно влияют на качественную картину развития процесса. 


При соединении металлических и неметалличе- 
ских материалов с использованием синтеза в твер- 
дой фазе, например, в процессе термитной сварки, 
возникает проблема выбора условий, обеспечи- 
вающих нужный режим синтеза. Соединение мате- 
риалов можно осуществить в разных режимах, за- 
висящих от технологических условий. При матема- 
тическом моделировании технологического про- 
цесса условия эксперимента определяют симме- 
трию и краевые условия. В соответствии с началь- 
ными и граничными условиями можно сформули- 
ровать различные задачи о соединении материалов: 
задачу о тепловом взрыве в толстостенном сосуде 
[1] и задачу о распространении фронта экзотерми- 
ческой реакции в узком слое между двумя инерт- 
ными материалами [2]. 

Подробное численное исследование [1-4] по- 
казало, что при соединении материалов с исполь- 
зованием синтеза в твердой фазе реализуются раз- 
личные режимы превращения. В задаче о тепловом 
взрыве в толстостенном сосуде численно исследо- 
ваны характеристики режимов превращения в за- 
висимости от параметров модели. В общем случае 
критические условия и время воспламенения зави- 
сят от размеров соединяемых материалов и толщи- 
ны слоя реакционной смеси, от теплофизических 
параметров инертного материала и реагента, а так- 
же от формально - кинетических параметров сум- 
марной химической реакции. Управление темпера- 
турой синтеза возможно за счет подбора соедини- 
тельного состава, доли и типа инертного наполни- 
теля в нем. Добавление в состав смеси инертного 
наполнителя позволяет реализовывать режим син- 
теза при заданной температуре. Дополнительная 
возможность снижения температуры синтеза по- 
является при использовании инертного наполни- 
теля с низкой температурой плавления. Роль плав- 
ления инертного наполнителя в соединительной 
смеси может быть существенной: от задержки вос- 
пламенения до смены механизма его влияния на 
воспламенение. Обнаружено, что, как без учета 
плавления, так и с ним, существует область пара- 
метров модели, где синтез осуществляется практи- 
чески при заданной температуре. Время превраще- 
ния в таком режиме на 1-3 порядка больше, чем в 
режиме теплового взрыва. 


С использованием метода малых возмущений в 
[5] показано, что доля инертного наполнителя не 
влияет на устойчивость стационарных режимов 
превращения при соединении материалов. Но, 
варьируя теплофизические свойства наполнителя 
(что в реальной ситуации соответствует замене од- 
ного наполнителя другим), можно добиться сниже- 
ния температуры синтеза соединительного шва, 
требуемого практикой. 

При численном исследовании задачи о распро- 
странении фронта экзотермической реакции в уз- 
ком слое между двумя инертными пластинами по- 
казано, что реализуются различные стационарные 
и нестационарные режимы превращения. В этом 
случае режимы превращения можно регулировать 
за счет варьирования теплофизических свойств ма- 
териалов. Численное исследование задачи с им- 
пульсным внешним источником показало [4] , что в 
такой системе наблюдаются самоподдерживаю - 
щиеся режимы превращения. Скорость распро- 
странения фронта зависит от соотношения тепло- 
физических свойств разных материалов и от пара- 
метров, характеризующих реакцию. Обнаружено, 
что самоподдерживающиеся режимы реализуются 
в более широкой области параметров, чем это сле- 
дует из классических представлений, что связано с 
влиянием соединяемых материалов. 

Особую роль в макрокинетике химических ре- 
акций, протекающих в твердой фазе, играют на- 
пряжения и деформации. Высокие температуры, 
достигаемые в зоне реакции часто нежелательны, 
поскольку сопровождаются большими термиче- 
скими напряжениями, приводящими к нежела- 
тельным последствиям и неустойчивости режимов 
превращения. Поэтому при моделировании техно- 
логического процесса дополнительно к изучению 
режимов протекания реакций, следует проводить 
оценку напряжений в ходе превращения, что и яв- 
ляется целью настоящей работы. 

1. Общая постановка задачи 

Напряжения, сопровождающие реакцию в про- 
гретом слое, могут быть различных типов: термиче- 
ские, концентрационные, (диффузионные, хими- 
ческие, фазовые, структурные). Термические на- 
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пряжения возникают в зоне реакции вследствие 
наличия высоких градиентов температур. Концен- 
трационные напряжения в общем случае связаны с 
реальным распределением концентраций реаген- 
тов и продуктов и меняются в результате физико- 
химических превращений. Структурные напряже- 
ния часто связывают с изменением пористости или 
с появлением системы трещин. В химически реаги- 
рующей системе фазовые, химические и структур- 
ные изменения связаны непосредственно с изме- 
нением температуры. Поэтому термическим на- 
пряжениям отводится наибольшее внимание. 

Рассмотрим модель возбуждения химической 
реакции, протекающей в слое между двумя инерт- 
ными материалами, с учетом напряжений, возни- 
кающих в прогретом слое. Пусть в узком слое меж- 
ду двумя инертными пластинами конечной толщи- 
ны помещено реагирующее вещество в смеси с 
инертным наполнителем, способное к экзотерми- 
ческому превращению. Полагаем, что механиче- 
ские свойства соединительного состава и инертных 
стенок близки (что и требуется получить в резуль- 
тате синтеза). Учет инертного наполнителя в сое- 
динительной смеси представляет отдельный инте- 
рес и является специальной задачей. 



Рис. 1. Иллюстрация к общей постановке задачи 



Рис. 2. К постановке задачи о тепловом взрыве: 1) реагирую- 
щая смесь, 2) инертные стенки 


В общем случае система уравнений, описываю- 
щих задачу о возбуждении реакции в слое между 
инертными материалами, включает уравнение те- 
плопроводности для реагирующего вещества: 


дТ, 

с, А д( -Л 
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уравнение теплопроводности для инертных мате- 
риалов: 
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уравнение химической кинетики для реакции, ко- 
торую можно описать схемой: 

А^В, (3) 

где А - суммарный реагент, В - суммарный про- 
дукт реакции: 

(4, 


и уравнения равновесия: 


гт , / ' °- 


(5) 


где с„ р„ Я,- - теплоемкости, плотности и теплопро- 
водности веществ, 7] - температура, (индекс 1 от- 
носится к реагирующему веществу, 2 - к инертным 
материалам), х— пространственная координата, і— 
время, а - степень превращения, Е а , (2 0 , к 0 - фор- 
мально-кинетические параметры, К - универсаль- 
ная газовая постоянная, от - компоненты тензора 
напряжений, связанные с деформациями обоб- 
щенными соотношениями 

с/оѵ. =2 /ц с1е ] + 8 у [Я (іе кк - КсПѴ\, (6) 
„ , 2 

К = Я + — ц - изотермический модуль всесторон- 
него сжатия, Я, ц — коэффициенты Ламе, % - пер- 
вый инвариант тензора деформаций, 8 Ѵ - символ 
Кронекера. Применимость (6) к сложным средам, 
содержащим внутренние поверхности раздела, 
обосновывается в работах [6, 1\. 

Если учитываются только термические напря- 
жения, то функция в соотношении (6) имеет вид: 
сПѴ = 3 а т ВТ. 

Если учитываются необратимые деформации, 
связанные с изменением свойств в ходе реакции, 
то функция имеет вид: 

сПУ = 3 [а т ВТ + (а в - а Л )Ва\, 


где а т - линейный коэффициент теплового расши- 
рения, а в , а А - коэффициенты концентрационно- 
го расширения суммарного продукта и реагента, 
определяемые в термодинамике аналогично коэф- 
фициенту теплового расширения: 


а = -!- 
" ЗК 


дѴ 
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В 2 т, С л .е 


где V— удельный объем тела, содержащего продук- 
ты с концентрацией С в =у и реагенты с концентра- 
цией С А = 1 —у. 

В литературе существуют оценки для коэффи- 
циентов концентрационного расширения для ин- 
дивидуальных веществ [8]; их можно определять и 
экспериментально теми же методами, что и а т . В 
случае упрощенного описания химического пре- 
вращения (3) разность (а в -а А ) отражает интеграль- 
ный характер реакции: если (а в -а А )> 0, то реакция 
идет с расширением объема; в противоположной 
ситуации - объем уменьшается. Для реакций тер- 
мического разложения, протекающих в прогретом 
слое в процессе зажигания, возможны обе ситуа- 
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ции. Для оценки разности коэффициентов можно 
использовать зависимость 


сс в а А 


= 1 Щ Ел__і 

3 [ р в т л 


где т,/ р,, і=А,В - молекулярные объемы реагента и 
продукта в своих фазах. 

В (5) не учитываются силы инерции, поскольку 
скорость распространения тепловых возмущений 
на 3-5 порядков меньше, чем скорость распростра- 
нения звука, и за время, много меньшее времени 
наблюдения (характерное для нормальных реакций 
в твердой фазе), волна механических возмущений 
успевает покидать зону прогрева. Эффектами свя- 
занности тепловых и механических процессов в 
уравнениях теплопроводности мы также пренебре- 
гаем, что оправдано оценками соответствующих 
коэффициентов и целью исследования. 

Математическая постановка задачи конкрети- 
зируется заданием условий возбуждения экзотер- 
мического превращения. 


2. Задача о тепловом воспламенении 
в сосуде с толстыми стенками 

Если в начальный момент времени соедини- 
тельная смесь и инертные материалы прогреты до 
одинаковой температуры Т п , причем Т 0 »Т е , где 
Т е - температура окружающей среды, то придем к 
задаче о воспламенении реакционной смеси в тол- 
стостенном сосуде. Примем, что размеры образца в 
направлениях осей О 2 и Ох много больше толщины 
слоя соединительной смеси, что позволяет прене- 
бречь теплопотерями в этих направлениях. В этом 
случае задача теплопроводности с химическими ис- 
точниками становится одномерной (рис. 2) и вклю- 
чает одномерное уравнение теплопроводности (1) 
(с производной по координате ) и уравнение хими- 
ческой кинетики (4) с начальными условиями: 

( = 0: Т г =Т 2 = Т 0 , а = О 


и условиями на границах: 


1 ду 2 ду ’ 1 - 


ВТ 

у = г + к:-Х 2 -ф = а(Т 2 - Т е ), 


где а - коэффициент теплоотдачи. 

Пусть торцы рассматриваемого толстостенного 
сосуда конечных размеров жестко закреплены по оси 
у, а остальные поверхности свободны от действия 
внешних сил, тогда в силу симметрии задачи имеем: 

Вуу = 0; е хх ф 0; ф 0 (7) 


< т а = 2 І* е „ + Хе кк -ЗКТѴ = 0 
■<т уу =Хе Ис -ЗШФ0 
оу = 2це ш + Хе № - 3 Ш = 0 


Данный случай соответствует плоской дефор- 
мации. Задача сводится к нахождению а я . 
Учитывая, что 8 кк =е хх +е а и е„=е а =е, найдем: 


8 


3 Ш 

2 р + X 


" °»~ гк№ -рк- 


( 8 ) 


Решение задачи проводится в безразмерных пе- 
ременных 


<7 

8 = ^-. 

а. 


Т -Т 

Ѳ= ±Е а 

КТ 2 а 


где <т* = 3 Ка т 



- характерное напряжение 


(максимальная величина напряжения, которая на- 
блюдается при нагреве до температуры Г 0 ). 

Уравнение (8) для искомого напряжения при- 
мет следующий вид: 

8 = -у т (Ѳ + Ѳ 0 + 8 ВА у), (9) 


где Ѳ 0 = — — т^Е - число Зельдовича, у т = — 

КТ 0 2 н + л 

а В4 Е 

Параметр § ВА = — - — - характеризует отношение 

(X у Р / ц 

химической деформации к термической. Если 
^>0, то реакция идет с расширением объема; в 
противоположной ситуации - объем уменьшается. 
Для реакций термического разложения, протекаю- 
щих в прогретом слое, возможны обе ситуации. 
Зная распределение температуры и концентрации, 
можем построить распределение напряжений. 

Задача о воспламенении решена численно по 
неявной разностной схеме с использованием мето- 
да прогонки. В расчетах определяются: время вос- 
пламенения, распределение температуры и степе- 
ни превращения, напряжения в различные момен- 
ты времени при варьировании параметров задачи, 
а также критические условия, разделяющие разные 
режимы превращения или «воспламенение» и «по- 
тухание» (т. е., прекращение реакции [1]). 

В качестве примера на рис. 2, 3 представлены 
пространственные распределения температуры (а) 
и термических напряжений (б) в последовательные 
моменты времени г=0,1; 0,25; 0,5; 0,75; 1; 1,1; 1,25 
для набора параметров і5Г с =^=1;Ѳ 0 =20;у=0,01 и 
различных значений параметра 



где ( а 


с, А ЕТ; 

тк 


ехр 



период адиабатической 


индукции; Х\ = ; г — радиус сосуда. Рис. 3 соот- 

с іРі 

ветствует «воспламенению», а рис. 4 - потуханию. 
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Легко видеть, что поведение напряжений пол- 
ностью определяется поведением температуры, 
что, разумеется, следует из (9) для § ВА =0. 




Рис. 3. Пространственное распределение температуры (а) и 
термических напряжений (6) в моменты времени 
г=0,1; 0,25; 0,5; 0,75; 1; 1,1; 1,25 для параметров 
К С =К Х = 1; 0 о =2О; у=0,01; <5=5; с=0; <5,=1,76 




Рис. 4. Пространственное распределение температуры (а) и 
термических напряжений (6) в моменты времени 
і=0,1; 0,25; 0,5; 0,75; 1; 1,1; 1,25 для параметров 
К С =К Х = 1; Ѳ 0 =20; у=0,01; 5=1; с=0; <5,=1,76 

Заметим, что в случае «потухания» степень пре- 
вращения ничтожна, а в случае воспламенения зо- 
на реакции на начальной стадии воспламенения 


достаточно узкая, чтобы концентрационные на- 
пряжения могли сказаться на величине напряже- 
ний существенным образом. Исключение соста- 
вляет область параметров модели, когда высокая 
теплопроводность соединяемых материалов при- 
водит к увеличению зоны реакции при воспламе- 
нении, или высокая их теплоемкость приводит к 
увеличению затрат тепла на их прогрев и увеличе- 
нию времени воспламенения, определяемому по 
резкому росту температуры в центре реагента. 

Иная ситуация складывается при добавлении в 
соединительную смесь инертного наполнителя. В 
математической формулировке задачи это приво- 
дит к введению понятия эффективных теплофизи- 
ческих свойств 

Кд =К( 1 - г) + Яг; с еЖ р еЖ = с 1 р 1 (1 - г) + ср г, 


где Хф Сф Рф - эффективные теплопроводность, 
теплоемкость и плотность веществ, % - доля инерт- 
ного наполнителя (индекс 1 относится к реагирую- 
щему веществу, величины без индекса характеризу- 
ют свойства инертного наполнителя в соедини- 
тельной смеси). 

На рис. 4, в, представлены пространственные 
распределения напряжений в разные моменты вре- 
мени с учетом инертного наполнителя ?=0,7 в сое- 
динительной смеси для различных значений § ВА для 
случая «воспламенения». Принято, что инертный 
наполнитель имеет те же теплофизические свой- 
ства, что и соединяемые материалы: 
К=В = \-КЛ=\. где 


К = 


с іРт 

с, А 




Я . п 

-ѵ 2 *- 


ср_ 

С,Р, 


■л - 


Х_ 

К' 


Добавление инертного наполнителя приводит к 
увеличению реакционной зоны и качественному 
изменению формы температурных кривых при 
воспламенении (рис. 2, а, и 4, а). Если § ВА =0, фор- 
ма кривых 5(<^,т) практически повторяет форму 
температурных кривых (рис. 4, б). Если реакция 
идет с увеличением объема (&Д0), напряжения в 
зоне реакции возрастают по абсолютной величине 
(рис. 4, в), что нежелательно для практического ис- 
пользования. Если реакция идет с уменьшением 
объема (§^<0), напряжения в зоне реакции суще- 
ственно меньше, чем в двух предыдущих случаях, 
но весьма велик перепад напряжений в зоне кон- 
такта соединительного состава и соединительной 
смеси (рис. 4, г). Разумеется, что в реальной техно- 
логической ситуации должен сформироваться не- 
который переходный слой между материалами в 
процессе синтеза или медленного остывания сое- 
динения, что в данной модели не исследуется. 

Зная свойства соединяемых материалов, можно 
провести соответствующие оценки напряжений. 
Так, если использовать в качестве соединительной 
смеси - смеси металлов, то в соответствии с лите- 
ратурными данными [9] имеем: 

Я=45...420 Па, а г =4,4...26 К 1 , 

ЕДТ20 кДжДмоль-К), Д=8,31 Дж/моль, Т 0 » 600 К. 
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а б 




в г 

Рис. 4. Пространственное распределение температуры (а) и напряжений (6-г) для набора параметров Ѳ 0 =5; у=5; с=0,03; <5=9; в 
разные моменты времени т=0,25 1; 2,5; 3,5; 5; 5,8; 6,4; 6,9; 7,1; а, б)$ 1А =0; в) % м = 5; г) ^ ВЛ =- 5 


Тогда, применяя известные соотношения: 

К = Е ■ .7 = ГУ _ Е 

3(1 — 2ѵ) ’ (1 — 2ѵ)(1 +ѵ) ’ ” 2(1 +ѵ) ’ 

можно найти параметры, входящие в модель. Пре- 
дел прочности оценим как оу^О, 1 Т упр . Тогда полу- 
чим область изменения максимальной величины 
напряжений при нагреве до температуры Г 0 (т. е. 
масштаб): 

сг,«14,8...816 МПа (10) 

и предел хрупкого разрушения или теоретический 
предел прочности: а пр ~4,5..А2 ГПа. 

Используя, характерный масштаб изменения 
напряжений (10), можем оценить величину макси- 
мальных напряжений в реакционной зоне, дости- 
гаемую в разных режимах синтеза (табл. 1). Такая 
оценка не учитывает многих факторов и является 
завышенной, но, тем не менее, позволяет сделать 
некоторые выводы. Так, в случае теплового взрыва, 
как с учетом инертного наполнителя, так и без на- 
полнителя в соединительной смеси напряжения в 
зоне реакции достаточно высокие. Для низкотем- 
пературного режима медленного превращения (ре- 


жима синтеза практически при заданной темпера- 
туре) [2], обнаруженного в модели с учетом инерт- 
ного наполнителя, напряжения в зоне реакции су- 
щественно меньше, что позволяет использовать 
этот режим для технической реализации процесса. 
Минимальное значение величины напряжений на- 
блюдается в некоторой области параметров для 
случая с инертным наполнителем в соединитель- 
ной смеси для реакции, идущей с уменьшением 
объема. В случае «потухания» напряжения также 
невысокие и не превышают теоретического преде- 
ла прочности. 

Очевидно, что существенную роль в величине 
напряжений в зоне реакции и на границе раздела 
материалов играют соотношения теплофизических 
свойств веществ. 

3. Оценка напряжений в задаче горения 

Если поджигание реагирующей смеси осущест- 
вляется с левого торца, где в момент времени 1=0 
температура поверхности х=0 повышается от Г 0 до 
Т,»Т 0 , то придем к задаче о распространении 


Таблица 1. Значения напряжений 5 для реакций, идущих с увеличением и уменьшением объема 


Режим превращения 

Без учета наполнителя 

С учетом наполнителя | 

1 

2 

3 

4 

5 

6 

йм=0 

&4> 0 

&м<0 

ІВА= 0 

&4>0 

йм<0 

Тепловой взрыв 

779.. .43000 

819.. .45206 

731. ..40310 

1181. ..65280 

1228. ..67728 

1139. ..62832 

Потухание (1-3) или медленное превращение (4—6) 

58. ..3223 

58. ..3223 

58. ..3223 

55. ..3035 

102. ..5630 

52.. .2856 
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фронта экзотермической реакции в узком слое 
между двумя инертными пластинами. Торцы рас- 
сматриваемой системы жестко закреплены по оси у 
(рис. 5), а остальные поверхности свободны от дей- 
ствия внешних сил, тогда соотношения для компо- 
нент тензоров напряжений и деформаций имеют 
вид, аналогичный (6). Т. е., опять имеем случай 
плоской деформации. Решение задачи о механиче- 
ском равновесии при условии постоянства свойств 
имеет вид (7, 8). Но в этом случае, в отличие, от 
предыдущего, тепловая задача становится двумер- 
ной и включает двумерное уравнение теплопровод- 
ности (с частными производными по осям Ох и Оу) 
с химическим источником и уравнение химиче- 
ской кинетики (4) С граничными условиями: 

дТ дТ 

У 0: Л 0: у^к/2 + Н-.-Ъ-ф = а{Т с -Т 2 )- 


ВТ дТ 

у = Ы2:^ = Ьф ^ 


х = 0 : 


Т = Т 1<1 

± І 


дТ 


дТ ; х -» / : = 0, і = 1, 2. 

_ с = 0 I > I ох 

дх 


В начальный момент времени имеем: 1=0: Т;=Т 0 \ 
У=0- 

Эта задача также решена численно в безразмер- 
ных переменных с использованием схемы расще- 
пления по координатам. 



Рис. 5. Иллюстрации к постановке задачи горения: 1) соеди- 
нительный состав; 2) инертные материалы 

На рис. 6 представлены изолинии (а) и двумер- 
ные поля напряжений (б) в разные моменты времени 
для набора параметров К с =\; К,= 0,1; 0 О =5 ; г=0,3; 
у=0,01 с учетом инертного наполнителя в соедини- 
тельной смеси. Аналогичные распределения пред- 
ставлены на рис. 7 для параметров К с = 0,1; А^=10; 
0 О =5; 7=0,0 1 ; ^=0,3. Из представленных рисунков 
легко видеть, что теплофизические свойства матери- 
алов значительно влияют на качественную картину 
развития процесса. Так, например, при низкой те- 
плопроводности соединяемых материалов, А^=0,1, 
фронт реакции опережает прогрев соединяемых ма- 
териалов, тепло из зоны реакции отводится плохо. 
Максимальные напряжения, тепловыделение и сте- 



Рис. 6. Изолинии (а) и двумерные поля напряжений (б) в моменты времени т=53, 60 для параметров К с =1, Лі=0,1; Ѳ 0 =5; у=0,01; г=0,3 
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Рис. 7. Изолинии (а) и двумерные поля напряжений (6) в моменты времени для параметров ^=0,1; Лі=10; Ѳ„=5; >=0,01; с=0,3 


пень превращения находятся в окрестности границы 
раздела веществ. В этом случае происходит посте- 
пенное (медленное) превращение вещества, кото- 
рое, тем не менее, поддерживается теплом, запасен- 
ным в инертных стенках. В иной ситуации вслед- 
ствие большой теплопроводности ^=10 и низкой 
теплоемкости К с = 0,1, инертные материалы прогре- 
ваются быстро, способствуя прогреву реагента: на 
инертной стадии прогрев соединяемых материалов 
идет быстрее. Тепла, запасенного в инертных стенках 
мало, но достаточно для протекания реакции. Вна- 
чале реакция идет по всей ширине соединительного 
слоя практически равномерно, распространяется на 
всю область постепенно, без ускорения и резкого те- 
пловыделения (чему способствуют толстые стенки), 
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затем реакция ускоряется и после длительной за- 
держки происходит тепловой взрыв, сопровождаю- 
щийся высокими напряжениями. 

Таким образом, в работе проведены оценки на- 
пряжений в зоне реакции, которые могут возни- 
кать в процессе соединения материалов с исполь- 
зованием синтеза в твердой фазе. Показано, что 
напряжения и деформации в соединительном со- 
ставе есть следствие высоких градиентов темпера- 
тур в зоне реакции и необратимого изменения 
свойств. Их величина в рамках макрокинетическо- 
го подхода полностью определяется распределени- 
ем температуры и степени превращения, как в ре- 
жиме горения, так и в режиме теплового взрыва. 
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